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Abstract 

Background: Shigatoxigenic E. coli are a global and emerging health concern. Shiga toxin, Stx, is encoded on the 
genome of temperate, lambdoid Stx phages. Genes essential for phage maintenance and replication are encoded 
on approximately 50% of the genome, while most of the remaining genes are of unknown function nor is it 
known if these annotated hypothetical genes are even expressed. It is hypothesized that many of the latter have 
been maintained due to positive selection pressure, and that some, expressed in the lysogen host, have a role in 
pathogenicity. This study used Change Mediated Antigen Technology (CMAT)™ and 2D-PAGE, in combination with 
RT-qPCR, to identify Stx phage genes that are expressed in £ coli during the lysogenic cycle. 

Results: Lysogen cultures propagated for 5-6 hours produced a high cell density with a low proportion of 
spontaneous prophage induction events. The expression of 26 phage genes was detected in these cultures by 
differential 2D-PAGE of expressed proteins and CMAT. Detailed analyses of 10 of these genes revealed that three 
were unequivocally expressed in the lysogen, two expressed from a known lysogenic cycle promoter and one 
uncoupled from the phage regulatory network. 

Conclusion: Propagation of a lysogen culture in which no cells at all are undergoing spontaneous lysis is 
impossible. To overcome this, RT-qPCR was used to determine gene expression profiles associated with the growth 
phase of lysogens. This enabled the definitive identification of three lambdoid Stx phage genes that are expressed 
in the lysogen and seven that are expressed during lysis. Conservation of these genes in this phage genome, and 
other Stx phages where they have been identified as present, indicates their importance in the phage/lysogen life 
cycle, with possible implications for the biology and pathogenicity of the bacterial host. 



Background 

Shigatoxigenic Escherichia coli (STEC) cause disease in 
humans following colonisation of the intestinal tract [1]. 
These infections are often serious, presenting with 
severe diarrhoea accompanied by haemorrhagic colitis. 
Downstream sequelae such as haemolytic uraemic syn- 
drome (HUS) and thrombotic thrombocytopenic pur- 
pura (TTP) can be fatal [2,3]. 

The principle defining virulence determinant of all 
STEC strains is the production of Shiga toxin (Stx), also 
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known as verocytotoxin (VT) or Shiga-like toxin (SLT) 
(1), of which there are two distinct forms, Stxl and Stx2 
[4]. Two variants of Stxl have been identified [5,6], 
whilst Stx2 is heterogeneous, with some variants more 
frequently associated with serious STEC outbreaks [1,7]. 
The stx genes are carried by temperate lambdoid bacter- 
iophages, which enter either the lytic or the lysogenic 
pathways upon infection of a bacterial cell [8-10]. Any 
bacteriophage encoding Stx is termed an Stx phage, and 
there is much genotypic and phenotypic diversity within 
this loosely-defined group [11]. Integrated Stx phages 
may exist in the bacterial chromosome as inducible pro- 
phages, or their residence within a host cell may facili- 
tate recombination events leading to the loss of 
prophage sequences, resulting in uninducible, remnant 
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Stx prophages within the lysogen chromosome [12]. The 
stx genes are located with genes involved in the lytic 
cycle; hence Shiga toxin expression occurs when Stx 
phages are induced into this pathway [11,13]. 

Stx phages possess genomes that are generally ~50% 
larger than that of the first described lambdoid phage, X 
itself, and -74% of Stx phage genes have not been defi- 
nitively assigned a function [11]. Genes that are essential 
for the Stx phage lifestyle are carried on approximately 
30 kb of DNA [14], whilst the entire genome is ca 60 kb 
in size in most cases [11,15,16]. The impact of Stx 
prophage carriage on the pathogenicity profile or biol- 
ogy of the host, beyond conferring the ability to produce 
Shiga toxin, has remained largely unexplored and it can 
be suggested that the accessory genome of Stx phages is 
likely to encode functions for which there has been 
positive selection [11]. 

In this paper, we describe the use of proteomic-based 
protein profile comparisons and Change Mediated Anti- 
gen Technology™ (CMAT) (Oragenics Inc.) [17] to 
identify Stx phage genes that are expressed during the 
lysogenic pathway. An E. coli lysogen of <t>24 B ::Kan, in 
which a kanamycin-resistance cassette interrupts the 
stX2A gene [18] of a phage isolated from an E. coli 
0157:H7 disease outbreak strain, was subjected to both 
CMAT and two dimensional polyacrylamide gel electro- 
phoresis (2D-PAGE) analyses of the expressed pro- 
teome. The <1>24 B ::Kan genome is 57.6 kb in size and is 
identical in all aspects to its wild-type parental phage 
other than the stxA gene interruption [14,18]. The 
majority of genes and coding sequences (CDS) carried 
by <D24 B are simply annotated as hypothetical [GenBank: 
HM 208303]. Bacteriophages tightly regulate expression 
of their genes involved in maintenance of lysogeny ver- 
sus replication of viral progeny, and the differentiation 
of gene expression associated with each state needed to 
be carefully determined in order to definitively associate 
expressed proteins and their genes with either the tem- 
perate or the lytic cycle. 

Results 

The rate of spontaneous lysis in an E. coli MCI 061 (<I>24 B ) 
culture at different stages of growth 

Spontaneous induction, defined as the induction of pro- 
phages from lysogens in the absence of an applied sti- 
mulus [19], occurs constantly in a proportion of the 
lysogen population in any culture, and this could ser- 
iously interfere with the differentiation of gene expres- 
sion between lytic and lysogenic states. In this study, it 
was necessary to determine culture conditions under 
which the number of spontaneous induction events was 
low whilst the cell density was high, enabling the consis- 
tent harvesting of sufficient amounts of cell-associated 
protein for downstream analyses. Lysogen cultures were 



sampled at hourly intervals beginning two hours post 
inoculation, and the c.f.u. ml" and p.f.u. ml" deter- 
mined. The lowest ratio of infective phages to cells, 
1:50, occurred at both 2 h and 3 h of lysogen growth. 
However the c.f.u. ml" during these times was relatively 
low; OD 600 = 0.184 (± 0.003) and OD 600 = 0.651 (± 
0.008), respectively. The ratio of phage to host cells 
increased sharply after 4 h of growth, before dropping 
after 5 h to 1:33 (OD600 = 1.192 [± 0.011]). The ratio 
of phage to cells in the culture remained stable at 1:33 
through to 6 hours of growth. Lysogen growth condi- 
tions were therefore standardised for MCI 061 (<t>24 B ) at 
5-6 hours when the cells were grown to an OD 60 o of 
1.2-1.3. 

Phage-encoded, lysogen-culture gene expression 
identified by CMAT 

A total of 13,519 clones were subjected to CMAT pri- 
mary screening, and taking efficiency of the library into 
account, this equates to a 3.3x coverage of the phage 
genome. Of these, 330 were identified by the lysogen- 
specific antiserum and chosen for further analyses and 
secondary screening. After two rounds of secondary 
screening, 250 clones were removed from the study and 
PCR analysis of the remaining 80 clones demonstrated 
that 46 possessed vector DNA only. The remaining 34 
recombinant transformants produced a peptide recog- 
nised by antibodies in the lysogen specific antiserum. 
The cloned inserts were sequenced, and the DNA 
sequences translated in all six possible reading frames. 
Twenty-three of the clones possessed sequences from 
twenty different <D24 B CDS (Table 1, Figure 1). The 
remaining eleven clones did not align with any <524 B 
-encoded CDS, although six did possess non-coding 
regions of the phage genome. The other five clones con- 
tained plasmid DNA only. 

Phage-encoded, lysogen-culture gene expression 
identified by 2D- PAGE 

Reproducible sets of gels from 2D-PAGE analyses were 
obtained through the utilisation of IPG strips in the pH 
ranges of 3.5-5.6 and 5.3-6.5. The optimal protein con- 
centration loaded on the gels was found to be 200 ug of 
total cellular protein from crude cell lysates. A total of 
42 protein spots were found only in the lysogen gel sets 
(data not shown); these were excised from the gels and 
analysed by MALDI-TOF. Twenty-four of these spots 
(Figure 2) contained enough protein for the generation 
of mass spectral data. When these spectra were searched 
against the University of Liverpool MASCOT database, 
which included all of the <1>24 B genome predicted pro- 
teins, six samples matched predicted phage proteins (PI 
to P6, Table 2, Figure 1). The remaining 20 spots were 
identified as E. coli proteins (Table 2); these are 
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Table 1 CDS identified by CMAT and location on the 
<t>24 B genome 



Clone 


Alignment to 024 B 
genome 


Aligned CDS 


Possible 
gene 


CM1 


39370-39772 


38090-40027 


tspS 


CM2 + 
CM14 


17489-18104 


17559-18086 


dam 


CM3 


2523-2185 


a: 2378-2286 
b: 2507-2379 




CM4 


3025-2375 


a: 2545-2375 
b: 2812-2711 
c: 2911-2840 




CM5 


54385-53866 


53693-53866 




CM6 


53690-53235 


53482-53297 




CM7 + 
CM13 


55160-55667 


49148-57571 




CM8 


38754-39248 


38460-38954 




CM9 


2542-2940 


2248-2646 




CM10 


35049-34598 


33695-34702 




CM11 + 
CM12 


39573-40016 


40189-39355 




CM15 


40137-40506 


40345-40626 




CM16 


38041-37623 


38000-37698 




CM17 


52465-52147 


52191-52514 




CM18 


45227-45877 


44818-45552 


lom 


CM19 


45610-46100 


45981-46382 




CM20 


4098-3676 


4333-4052 




CM21 


39305-39919 


39405-39650 




CM22 


39875-40526 


39909-40298 




CM23 


45713-46232 


a: 45784- 
45921 

b: 46072- 
46239 





potentially lysogen specific but were not investigated 
further here. 

Analyses of gene expression patterns 

Generally, lambdoid phage regulatory circuits tightly 
control the expression of genes, yet some of the genes 
identified in the CMAT library and the 2D-PAGE ana- 
lyses above were phage genes whose expression should 
be linked to prophage induction (Figure 1) and not the 
stable prophage state, e.g. the gene encoding the tail 
spike protein. It was assumed that gene expression nor- 
mally linked to the lytic replication cycle must be at a 
very high level in a small subset of the cells and that 
lysogen-restricted gene expression patterns of these 
genes might be very low, especially as neither CMAT 
nor 2D-PAGE identified the expression of repressor, the 
product of the cl gene, in the lysogen culture. Therefore 
it was essential to devise a method that would deter- 
mine whether phage genes were being expressed by the 
majority of the stable lysogen population, or the small 
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Figure 1 Schematic representation of the 024 B genome 

Squares symbolise the locations of the CMAT and PAGE CDS 
identified as well as some of the essential genes involved in the life 
cycle of the phage. - represents 5 kb. For further details on the 
gene identities see Tables 1 & 2. 
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4 pH 7 3 OH 5.6 

Figure 2 2D-PAGE images of total cell protein from MC1061/O24 B ::Kan IEF on pH range 4-7 (A, C), 5.3-6.5 (B) and 3-5.6 (D). Arrows 
represent proteins identified as phage encoded; circles represent proteins identified as encoded by £ coli, but not present on corresponding 
naive MC1061 gels (data not shown). 



subset of the population undergoing spontaneous induc- 
tion events. A strategy involving qPCR was developed to 
provide this important information, and a variety of 
genes were chosen as controls for this study (Additional 
file 1: Table SI, Figure 1). Calibration curves for quanti- 
tation and comparison of the qRT-PCR data were pro- 
duced for every set of primers used; R values from 
linear regression analyses of these standards ranged 
from 0.990 to 0.999 with slopes ranging between -3.72 
and -3.10 (Additional file 1: Table SI). 

The data from the qPCR assay were analysed by com- 
paring the shape of the expression data for any given 
gene from a lysogen culture throughout the prophage 
induction process where time 0 is the point of norfloxa- 
cin (inducer) addition (Figure 3). Lysogen-restricted 
gene expression should be negatively affected after 
induction (Figure 3A, CI), and if expression is actually 
linked to the small proportion of cells undergoing spon- 
taneous induction, then the expression levels should rise 
during the induction process. This is indeed the case as 



expected for Q, Cro, Capsid & Terminase, which display 
a significant increase after 50 min of recovery, Figure 3; 
Additional file 2: Table S2). 

Four genes identified by 2D-PAGE, PI, P4, P5 and P6, 
visibly follow the same expression pattern as the genes 
expressed during the lytic cycle and accordingly the 
increase in gene copy number is significant (p-value < 
0.05) after 50 or 60 min of recovery from exposure to 
norfloxacin (Figure 3; Additional file 2: Table S2). P3 
and P2 appear to have a similar pattern to cl, i.e. their 
levels of expression in the lysogen are higher than the 
levels after induction; however the ANOVA analysis did 
not identify these differences as significant, probably due 
to the high variability amongst the replicates. Of the five 
CDS identified by CMAT, which were subsequently 
selected for expression analysis based upon either their 
putative function or location within the phage genome, 
four had expression patterns linked to the lytic cycle. 
CM18 was shown by qPCR to be strongly expressed in 
lysogen cultures, but when the cells are induced, high 
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Table 2 Protein identities according to the MASCOT database 





Gene name 


Access No. 


pi/ 

MIA/ 

(Da) 


Description 


Sequence 3 
Coverage 

(%) 


MASCOT b 
Score 


Peptides' 1 
matches 


Estimated 
ni/n/iiA/ rr*n\ 


MASCOT 
Database 
Identified in 


1 


PI 




5.28/ 
33860 


Identical to hypothetical 
protein p78 from 933 W d 


32 


63* 


6 


5.50/40000 




2 


P2 




5.27/ 
17096 


Similar to hypothetical protein 
p23 from 933 W 


42 


39 


5 


r r\r\ /1 rr\r\r\ 

5.00/ 1 5000 




3 


P3 




5.09/ 
13472 


Similar to hypothetical protein 
p24 from 933 W 


33 


55 


3 


5.6/8000 




4 


P4 




5.14/ 
25800 


Identical to exo of 933 W 


29 


52 


5 


5.5/40000 




5 


P5 




5.29/ 
7336 


Not known homologue 


83 


■11 


2 


7.00/7000 




6 


P6 




5.22/ 
13628 


Identical to hypothetica 
protein Stx2lp064 e 


38 


30 


4 


7.00/1 0000 


1 


7 


nanA2 


Q6KD26 


5.77/ 
34077 


N-acetylneuraminate lyse 2 


21 


47 


'1 


5.3/35000 


2 


8 


gadB 


CAQ31981 


5.29/ 
52634 


Glutamate decarboxylase beta 


23 


57 


7 


5.3/35000 


2 


9 


sodB 


P0AGD5 


5.58/ 
21121 


Superoxide dismutase [Fe] 


■-10 


53 


6 


4.00/1 00000 


2 


10 


napA 


AAC75266 


8.23/ 
92983 


Nitrate reductase 


I 1 


49 


7 


5.5/100000 


2 


11 


tig 


AAA62791 


4.73/ 
47994 


Trigger factor 


24 


58* 


7 


3.5/6000 


2 


12 


UTI89_C3021 


Q1R837 


4.70/ 
6971 


Hypothetical protein 


70 


42 


2 


5.5/7000 


2 


13 


2FPKA 


ZP_04873224 


5.24/ 
32497 


6-phosphofructokinase 


23 


46 


5 


5.3/50000 


2 


14 


gcpE 


S23058 


5.87/ 
40658 


1 -hydroxy-2-methyl-2-(E)- 
butenyl 4-diphosphate 
synthase 


16 


38 


-1 


5.5/80000 


2 


15 


aceE 


P0AFG9 


5.46/ 
99475 


Pyruvate dehydrogenase E1 
component 


10 


60* 


7 


5.4/100000 


2 


16 


bfpK 


Q9S141 


7.63/ 
18294 


BfpK 


5-1 


49 


3 


6.4/25000 


2 


17 


ychN 


P0AB53 


5.02/ 
12685 


ychN 


'16 


38 


2 


5.3/100000 


2 


18 


UTI89_C1 147 


Q1 RDD 


5.57/ 
24945 


Hypothetical protein 


15 


38 


'1 


5.7/30000 


2 


19 


ompC 


Q9RH85 


4.55/ 
40474 


Outer membrane protein 


18 


'11 


>\ 


5.5/40000 


2 


20 


ECs1247 


G90784 


4.74/ 
25144 


Hypothetical protein 


26 


39 


5 


6.5/35000 


2 


21 


UTI89_C2748 


Q1R8V6 


10.19/ 
10724 


Hypothetical protein 


44 


50 


-1 


6.4/8000 


2 


22 


nirB 


E86001 


5.79/ 
93112 


Nitrate reductase (NAD(P)H) 
Subunit 


10 


53 


8 


5.3/1 00000 


2 


23 


yagP 


CAQ30761 


5.65/ 
15401 


yagP protein 


36 


43 


3 


5.3/10000 


2 


24 


rhsF 


Q47284 


5.69/ 
23342 


RhsF 


18 


42 


-I 


5.69/8000 


2 



Table represents matches to f. coll proteins in the MASCOT database and matches to 024 B proteins in the University of Liverpool local MASCOT database 
a percentage of sequence of the matched protein that is covered by the experimental MS. 

b logarithm of the probability that the match between the experimental data and a protein sequence in the database is a random event. 
c number of peptides that match the protein in the database 
d 933 W is an Stx2 phage described by Plunkett ef al. [16]. 
e Stx2 is an Stx2 phage described by Sato ef al. [20]. 
^represents significant matches (p-value < 0.05) 

1 University of Liverpool local MASCOT database; 2 general MASCOT database 
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Figure 3 Graph depicting gene expression profiles before and following norfloxacin induction. Panel A: Control Genes. CI, marker gene 
for lysogeny-restricted expression; Cro and Q, marker genes for early induction response; Ter and Cap, marker genes for late gene expression; 
GyrB and 16S, marker genes for the cellular response. Panel B: Expression profiles of the prophage genes identified by CMAT. Panel C: 
Expression profiles of the genes identified by 2D-PAGE analyses of the lysogen. The Y axis represents gene copy number per 300 ng of RNA; the 
X axis represents time (min). Time -60 refers to the samples taken before induction and represents the lysogen population, Time 0 represents 
samples taken at the beginning of the recovery time, Time 10, 10 min after recovery, etc. The experiment was run using biological replicates, but 
due to the asynchronicity of induction across these experiments the data from a representative single biological replicate are shown. 



expression levels are maintained, suggesting that expres- 
sion of this gene has been uncoupled from the phage 
regulatory circuits. The outcome of one-way ANOVA 
analysis to determine the impact of prophage induction 
on gene expression was found to be significant in 11 
cases (p-value < 0.05): cl, cro, terminase, capsid, Q, 
CM1, CM2, CM5, CM7, PI and P5. The other 7 genes 
studied did not present significant changes in expres- 
sion: P2, P3, P4, P6, CM18, 16S, and gyraseB. The full 
set of p-values for the data in Figure 3 are presented in 
Additional file 2: Table S2. 

Discussion 

Temperate phages, maintained as prophages in their 
lysogens, have been the subject of speculation concern- 
ing their benefit to the host: selective advantage, 
increased virulence, and other traits with varying 
degrees of direct and/or indirect impact on the host 
have been identified [11,21-27]. The challenge in this 
area has been how to identify phage-encoded genes that 
directly affect their lysogen, because many/most phage 
genes are annotated as encoding hypothetical proteins. 
In addition, there will always be a small background 
population undergoing spontaneous induction in the 
absence of discernible stimuli [19], potentially confound- 
ing the identification of lysogen-restricted prophage 
gene expression. In a specific E. coli lysogen of Stx2- 
phage 933W, a phage very closely related to 024 B , the 
spontaneous induction rate was calculated as 0.014% 
[28], which means that in a lysogen culture fourteen 



cells per 100,000 are undergoing prophage induction. 
Other recent work was demonstrated that various 
induction agents and growth conditions differentially 
effects induction in a prophage-dependent manner [29]. 
Assuming a burst size similar to that of bacteriophage 
Lambda (170 ± 10 virions cell" 1 ) [27], a significant 
amount of phage structural protein production can 
occur in an uninduced lysogen culture. 

In order to mitigate this effect, the growth phase at 
which the ratio of lysogens to free phage was high (two 
to three hours post inoculation) was targeted. However, 
the cell density at this point was very low and 5-6 hours 
was chosen as the standardised incubation time as a 
compromise. In this study, 26 genes from the bacter- 
iophage 024 B were identified by either CMAT or 20- 
PAGE as being expressed in E. coli lysogen culture. No 
genes were identified by both CMAT and 2D-PAGE 
methods, perhaps due in part to the low absolute num- 
ber of <D24 B genes identified by the latter approach. 
However, the level of redundancy in the genes identified 
by the CMAT clones was lower than expected, given 
the number of clones screened and the calculated phage 
genome coverage; however, putative positive clones were 
selected conservatively in an attempt to limit the num- 
ber of false positives. Additionally, CMAT-based identi- 
fication may also introduce bias into library screening 
due to differences in protein immunogenicity and anti- 
genicity. It is important to note that the best charac- 
terised lysogen-restricted gene, cl (encoding lambdoid 
phage repressor), was not identified using either CMAT 
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or 2D-PAGE, indicating that this study was not exhaus- 
tive. Nevertheless, the paucity of information on lyso- 
gen-restricted gene expression is such that these data 
represent a significant step forward in our understand- 
ing of phage/host interactions and lysogen biology. 

Of the 26 phage genes identified in this study, Tsp, 
encoding the characterised tail spike protein of <1>24 B 
[30,31] was a known structural protein and therefore 
not expected to be expressed by a stable lysogen 
(Tables 1 &3), while the expression profiles of the 
other 25 proteins were unknown. Therefore the result- 
ing challenge was to identify the fraction of the culture 
(lysogens or cells undergoing lysis) that were responsi- 
ble for expression of these 26 phage genes as well as 
determining testable hypotheses to assign function to 
the identified gene products. Five genes identified dur- 
ing the CMAT screening were chosen for gene expres- 
sion profiling due to their genome location, potential 
function or degree of conservation across a range of 
phages (Table 3). The CDS CM18 encodes a Lorn 
orthologue, which was expected to be expressed in the 
lysogen as the lambda lorn gene is associated with the 
alteration of the lysogen's pathogenic profile after loca- 
tion of Lorn in the outer membrane [32-34]. However, 
expression of lorn in the <D24 B lysogen unexpectedly 
appears to be uncoupled from the phage regulatory 
pathways, because it is expressed at similar levels in an 
infected cell regardless of whether that cell exists as a 
stable lysogen or is undergoing prophage induction. 
The CDS CM2 encodes a putative Dam methyltransfer- 
ase. Bacterial-encoded Dam methyltransferase has been 
shown to be essential for maintenance of lysogeny in E. 
coli infected with Stx-phage 933 W [35]. The expres- 
sion pattern of the 024 B -encoded Dam methyltransfer- 
ase could indicate that it is fulfilling a similar role, or 
supplementing the function of the host-encoded Dam 
methylase in lysogens infected with this phage. The 
functions of CM5 and CM7 are unknown. CM7 is an 
ORF of 8 kb, and as the amount of DNA that can be 
packaged by a phage is limited, such a large gene is 
likely to be conserved only if it confers an advantage to 
the phage or its lysogen; it may be significant that this 
large gene is associated with several other phages 
(Table 3). CM5 is a small CDS located on the comple- 
mentary strand to the one encoding CM7, in a region 
with few other CDS, though it is directly upstream of 
another CMAT-identified CDS, CM6. The data (Figure 
3) indicate that the expression of these 3 genes is 
linked to prophage induction, a surprising outcome as 
CM7 does not appear to be a phage structural gene, 
has been indicated by bioinformatic analyses (data not 
shown) to be a probable outer membrane protein, and 
is downstream of CM18, whose regulation is uncoupled 
from expression of the late genes. 



The qPCR expression profile for the phage genes iden- 
tified as being expressed in the lysogen by 2D-PAGE, 
PI, P2, P3, P4, P5 and P6, indicated that only the 
expression of P2 and P3 were restricted to lysogen cul- 
tures with a stable prophage. The genes for both P2 and 
P3 lie downstream of the cl gene. However, their 
expression levels are one and five orders of magnitude 
greater, respectively, than the expression levels of cl, the 
lambdoid phage repressor gene. It is known that in 
Lambda phage, the cl gene transcript is leaderless, pos- 
sessing no ribosome binding site for initiation of transla- 
tion, with transcription and translation beginning at the 
AUG start codon [36]. If this causes the 5' end of the 
transcript to be less stable and more easily subject to 
degradation, the higher level of P3 transcript could sim- 
ply be due to possession of a longer half life than those 
genes at the 5' end of the transcript. 

The genes encoding P2 and P3 are conserved in many 
other phages (Table 3). They have no bioinformatically 
identifiable promoters of their own, so are likely to be 
driven by pRM or pRE like cl (see [37] for a cogent 
review of the related lambda phage), but differences in 
the levels of transcription between these 3 genes implies 
that there is still more to discover about the right opera- 
tor region of this phage. The proteins PI, P4, P5 and P6 
all exhibit gene expression profiles that suggest they are 
expressed following prophage induction. These genes 
are scattered across the phage genome (Figure 1) and 
are shared by various phages (Table 3). The protein P4 
appears to be part of the lambda Red recombinase sys- 
tem [38-40] and the data presented here suggest that 
this is most active upon prophage induction. This could 
be relevant to the mechanisms that underpin diversifica- 
tion, evolution and production of new phages by lyso- 
gens carrying an inducible prophage along with one or 
more inducible or remnant prophages [11,41,42]. The 
proteins PI, P5 and P6 are scattered across the genome 
on the strand typically associated with expression of 
genes linked to lysogenic infection (e.g. cIII, N, cl). Two 
genes encoding proteins PI, P5 and P6 are found in 
other phages, but have no known function. 

In summary, genome sequencing of prophages and 
bacteriophages has identified that these viral elements 
encode higher numbers of hypothetical genes than those 
to which we can currently assign a function. These 
genes are often conserved across many bacteriophages, 
but do not appear to encode structural proteins. For 
these genes to remain present in the phage genome, 
especially considering the fluidity of the genetic compo- 
sition of lambdoid phages [43], they must surely provide 
an important function in either the phage life cycle or 
that of the lysogen itself. In attempting to identify 
prophage genes whose expression was restricted to the 
stable prophage state, our goal was to identify prophage 
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Table 3 Distribution of the proteins identified by CMAT and 2D-PAGE across phage genomes 


fZ, Ann nthoi" Ctv nh^n^c rsrrwinn tno y\ rntfli nc in trio eti in\/ fifiantitwl 
Ucllc XJlllfcri jIa pildljtrb v-aliyllllj Lllc fJIUltrllii 111 lllc blUUy ^lUcllLlty/ 


r\(_Lfcrbil Ul 1 WllleT fJlldytri 


A rracc i/~in 
MLlcbblUll 




number 


number 


CM1 Stx2 converting phage II (99%) 


YP_003828920.1 




phage VT2-Sakai (99%) 


NP_050557.1 




phage 933 W (99%) 


NPJ349519.1 




Stx1 converting phage (99%) 


YP_003848832 




phage BP-933 W (99%) 


YP_003848832.1 




phage VT2phi_272 (99%) 


ADU03741.1 




phage Min27(100%) 


ADU03741 




CM2 Stx2 converting phage II (100%) 


BAC781 16 




phage VT2-Sakai (100%) 


NP 050531.1 




phage Min27 (100%) 


YP 001648926 




phage HK97 (99%) 


AAF31 1 37 




phage Lahn2 (99%) 


CAJ26400 




phage Lahn3 (98%) 


CAC95062.1 




phage 2851 (99%) 


CAQ82016 




phage CP-1 639(99%) 


CAC83142 




prophage CP-933 V(99%) 


AAG57233 




Phage Nil2 (99%)(99%) 


CAC95095 




Stx1 converting phage (99%) 


YP_003848889.1 




Phage CP-1 639 (99%) 


CAC83 142.1 




Phage YYZ-2088 (99%) 


YP_0022741 70.1 




Stx2-converting phage 1717 (99%) 


YP_002274244.1 




CM5 phage Min27 (100%) 


YP_00 1648966.1 




Stx2 converting phage 11(100%) 


YP_00388933.1 




Stx2 converting phage 1(100%) 


NP_6 12929.1 




phage VT2-Sakai (100%) 


NP_050570.1 




phage 933 W (1 00%) 


NP_049532.1 




phage VT2phi_272 (100%) 


ADU03756 




CM7 phage VT2-Sakai (99%) 


NP_050570 




Stxl converting phage (99%) 


BAC77866.1 




Phage VT2phi_272 (97%) 


ADU03756.1 




Phage 933 W (97%) 


NP_049532.1 




Stx2 converting phage 1 (97%) 


NP_6 12929.1 




Stx2 converting phage 11(97%) 


BAC78032.1 




Phage BP-933 W (97%) 


AAG55616.1 




Stx2 converting phage 86 (91%) 


YP_794082.1 




Phage Min27 (97%) 


YP_00 1648966.1 




CM18 phage VT2-Sakai (100%) 


NP_050564.1 




Stxl converting phage (100%) 


YP_003848839.1 




Phage 933 W (100%) 


NP_049526.1 




Stx2 converting phage 1 (100%) 


ZP_02785836.1 




Stx2 converting phage II (100%) 


YP_003828926.1 




Phage BP-933 W (100%) 


NP_286999.1 




Stx2 converting phage 86 (97%) 


YP_794076.1 




Phage Min27 (100%) 


YP_001 648959.1 




PI Stx2 converting phage II (99%) 


YP_003828937.1 Phage phiV10 (78%) 


YP_5 12303.1 


Stx2 converting phage 1 (99%) 


NP_61 2952.1 




Phage 933 W (99%) 


NP_049538.1 




Phage BP-933 W (99%) 


AAG55619.1 




phage VT2-Sakai (99%) 


NP_050575.1 




Phage Min27 (96%) 


YP_00 164890 1.1 
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Table 3 Distribution of the proteins identified by CMAT and 2D-PAGE across phage genomes (Continued) 





Stx2-converting phage 86 (96%) 
Phage BP-4795 (96%) 
phage CP-1639 (74%) 


YP_794094.1 
YP_00 1449244.1 
CAC83 133.1 




P2 


Stx2 converting phage I (100%) 
Phage 933 W (1 00%) 
Phage BP-933 W (100%) 
Phage Min27 (100%) 
Stx2-converting phage 86 (100%) 


NP_61 2997.1 
NP_049484.1 
AAG55573.1 
ABY49878.1 
YP_794 109.1 


Salmonella enteric YP_002455860.1 
bacteriophage SE1 (86%) 
Salmonella phage ST160 (86%) YP_0041 23782.1 


P3 


Stx2 converting phage 1 (100%) 

Phage 933 W (100%) 
Stx2-converting phage 86 (100%) 
Phage Min27 (100%) 
Phage BP-933 W (100%) 


NP_6 12995.1 
NP_049483.1 
YP_794 108.1 
YP_001 64891 5.1 
AAG55572.1 




P4 


Phage 933 W (1 00%) 

Phanp RP-Q33 W H OfWnl 
nidge Dr yjjj vv 1 1 uu/oj 

Prophage CP-933 V (100%) 

Stx2 converting phage 1 (100%) 

Phage VTI-Sakai (100%) 

Phage YYZ-2008 (99%) 

Stx2-converting phage 1717 (98%) 

prophage CP-933 K (98%) 

phage BP-4795 (98%) 

phage Min27 (99%) 


NP_049473.1 
MP ?R^Q^7 1 

NP_288695.1 
NP_6 12980.1 
BAB19617.1 
YP_002274150.1 
YP_002274221.1 
YP_003500773.1 
YP_00 1449249.1 
YP_001 648905.1 


Phage lambda (98%) NP_040616.1 


P5 


Stx2 converting phage I (100%) 

Phage 933 W (100%) 
Stx2 converting phage II (100%) 
Stx2-converting phage 1717 (98%) 
phage 2851 (98%) 
Phage BP-4795 (97%) 


NP_61 3032.1 
NP_049503.1 
BAC78139.1 

YP_002274255.1 
CAE53952.1 

YP_00 144 19282.1 




P6 


Stx2 converting phage 1 (99%) 
Stx2 converting phage II (99%) 
phage VT2phi_272 (99%) 

phage Min27 (99%) 
phage VT2-Sakai (99%) 
Stxl converting phage (99%) 
Stx2-converting phage 86 (96%) 


NP_61 2943.1 
BAC78046.1 

ADU03756.1 
YP_00 1648966.1 

NP_050570.1 
BAC77866.1 
BAF34067.1 





genes that were candidates for influencing the fitness of 
the bacterial host. However, the study was hampered by 
the fact that lysogen-restricted gene expression can be 
at very low levels, and phage genes associated with 
phage replication are expressed at very high levels. 

Conclusions 

Two different experimental strategies were employed to 
identify prophage genes expressed by their lysogen, and 
it is interesting to note that lysogen-specific antibody 
recognition of a peptide expression library and differen- 
tial 2D-PAGE with subsequent protein identification by 
peptide mass spectrometry, did not identify the same 
genes or proteins. The failure of both to identify 



expression of the cl gene encoding the phage repressor 
was shown by RT-qPCR to be due to the very low 
expression levels peculiar to this phage gene (Figure 4); 
the CI protein is also very susceptible to autocatalysis 
and therefore elusive. Both CMAT and 2D PAGE identi- 
fied some phage genes that were associated with lytic 
induction, and the qPCR strategy was useful for discri- 
minating low level expression in stable lysogens from 
high-level gene expression in the minority of lysogens 
that were undergoing spontaneous induction. Improving 
our understanding of the STEC disease process is ever 
more urgent in light of the recent emergence of a new 
Shiga-toxin producing E. coli pathotype [44], and deter- 
mining the function and expression patterns of the 
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genes in Stx phage genomes is very important in that 
context. 

Methods 

Bacterial strains and culture 

The E. coli K-12 strain, MC1061, was used as the bac- 
terial host for the production of lysogens. MC1061 
(024 B ) refers to the <1>24 B lysogen of MC1061; naive 
MC1061 refers to cells that have not been infected by 
<524 B . E. coli K-12 strain DM1187 was used as the indi- 
cator host strain in plaque assay experiments [18]. 
BL21-AI cells (Invitrogen, Paisley, U.K.) were used as 
the expression host for genetic constructs. Bacterial 
strains, plasmids and phages used in this study are listed 
in Table 4. 

All cultures, unless otherwise stated, were propagated 
from an overnight (-16 h) starter culture (0.5% v/v 
inoculum) in Luria Bertani (LB) broth (Merck KGaA, 
Darmstadt, Germany) containing 0.01 M CaCl2, incu- 
bated at 37°C with shaking at 200 r.p.m. Lysogen cul- 
tures were grown in the presence of kanamycin (Kan, 50 
ug ml" 1 ). Induction of protein expression in BL21-AI 
cells took place in BHI broth with 0.2% arabinose and 1 
mM IPTG. 

Induction of phage lysogens 

Cultures of MC1061(<X>24 B ) cells were incubated with 
norfloxacin (1 ug mL' ) for 1 h at 37°C with shaking at 
200 r.p.m. Cultures were then diluted 1:10 in fresh LB 
and the bacteria allowed to recover from the growth 
inhibitory effects of the antibiotic for 1 h at 37°C (the 
recovery period), with shaking at 200 r.p.m. 

Antisera production for use in CMAT 

A 2 L culture of MC1061(O24B) was propagated for 6 
hours. The cells were pelleted and resuspended in 1 ml 
of retained supernatant plus 1 ml of LB broth. Protease 
inhibitors (20 u.L) (Roche Complete Mini EDTA Free 
protease inhibitor cocktail tablets, Bath, U.K.) and 10 uL 
of lysis buffer (7 M urea, 2 M thiourea, 2% CHAPS, 1% 
DTT, Roche Complete Mini EDTA-free protease inhibi- 
tor cocktail tablets) were added to each. The samples 
were sonicated at 15-18 u for 6 x 10 s bursts. Absolute 
methanol (1.5 ml) was added, and the samples were 



incubated at -20°C for 60 min. Protein was harvested by 
centrifugation at 16,000 g for 5 min, and the resultant 
protein pellets were air-dried and suspended in 0.5 ml 
phosphate buffered saline (PBS). The samples were 
pooled; the protein content was measured by Bradford 
Assay [46] and adjusted to 1 mg ml" 1 . A total of 4 mg 
of the lysogen protein was sent to Eurogentec (Seraing, 
Belgium) for antisera production in rabbits, using the 
Ribi adjuvant system. Two rabbits were immunised with 
the protein sample on days 0, 14, 28 and 56 of the pro- 
gram. Bleeds were carried out on days 0 (pre-immune 
sera), 38, 66 and 87 (final bleed). Pre-immune sera from 
the two rabbits used were received and tested for cross- 
reactivity by western blot analysis. 

CMAT was carried out as per instructions from the 
license holder, Oragenics Inc., FL., U.S.A. [17,47], with 
the exception that BL21-AI was used as the expression 
strain for the phage library. The recommended expres- 
sion host, BL21[DE3], is an E. coli-'k lysogen, and there- 
fore an inappropriate strain to use in phage protein 
expression studies [48]. The expression library was cre- 
ated from <D24 B ::Kan DNA. The rabbit antisera were 
depleted of antibodies reactive to E. coli proteins by a 
series of adsorptions to naive MC1061 whole cells and 
cellular lysate, and to BL21-AI + pET30c (empty vector) 
whole cells and cellular lysate. The depleted antisera 
were compared to undepleted antisera by western blot. 
Adsorptions were repeated until no bands were detect- 
able by western blot probing of 6 (ig of naive MC1061 
proteins. 

Peptide expression library construction 

Semi-confluent plaque assay plates [18] were overlaid 
with 3 ml SM buffer (100 mM NaCl, 8 mM MgS04, 50 
mM Tris-HCl, pH 7.5) and incubated at 4°C for 16 h, 
with gentle agitation. The SM buffer and top agar were 
transferred to separate 50 ml centrifuge tubes that were 
vortexed with 10% (v/v) fresh SM buffer and subjected 
to centrifugation at 10,000 g for 10 min. The superna- 
tant was pooled and 30 (il of chloroform were added to 
each 10 ml of buffer. DNase (5 ug ml" 1 ) and RNase (1 
mg ml" ) were added, and the samples were incubated 
at 37°C for 1 h. PEG 8000 (33% [w/v]) was added, and 
the samples were incubated on ice for 30 min. 



Table 4 Bacterial strains, plasmid and phages used in the study 


E.coli strains, plasmids and phages 


Relevant Genotype 


Reference 


BL21-AI 


F- ompT hsdSB(tB-, mB-) gal dcm (DE3), arabinose inducible T7 RNA polymerase 


Invitrogen, Paisley, U.K. 


MC1061 


F" A(ara-/eu)7697 &{codB-lad)3 galK16 X' mcrAO rps/J50(strR) mcrBI 


[18] 


DM1187 


F- dam-13::Tn9(Cm R ) dcm- mcrB hsdR-M + gal! am- lac- thr- leu- tsxR 


[45] 


TOP10 


F- mcrA <D80/acZAM15 recA + 


Invitrogen, Paisley, U.K. 


pCR®-Blunt 


lacZ a, Kan R , ccdB 


Invitrogen, Paisley, U.K. 


pET30c 


Expression vector with T7 promoter, Kan R , Tet R , 


Novagen, Notts, UK 


<S24 B 


Stx2-phage, AstxA 2 "£iph3 


[14] 
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Precipitated phage particles were harvested by centrifu- 
gation for 10 min at 10,000 g, and the pellets were 
resuspended in 500 ul SM buffer per 30 ml starting 
volume. Samples were treated with DNase and RNase, 
as before. Phage DNA was purified by phenolxhloro- 
form:isoamyl alcohol extraction and isopropanol precipi- 
tation [49] and resuspended in 100 ul ddH 2 0. The 
024 B DNA (15 ug ml" in TE) was fragmented using a 
HydroShear (GeneMachines, MI, USA), at speed code 6 
for 30 cycles, followed by 30 cycles at speed code 2. 
DNA of the required size range (300-900 bp) was iso- 
lated by gel purification. pET30c plasmid (EMD Bios- 
ciences) DNA was digested with EcoR V and 
dephosphorylated with calf intestinal phosphatase (New 
England Biolabs) according to the manufacturer's 
recommendations. The size fractionated <I>24 B DNA 
fragments were cloned into the prepared pET30c DNA 
(50 ng) vector in a molar ratio of 25:1 (insert to vector). 
Chemically competent BL21-AI expression host cells 
(Invitrogen) were transformed with the plasmid DNA 
according to the manufacturer's recommendations. 

Primary screening 

Transformed BL21-AI cells were plated onto LBKan 
plates and incubated at 37°C (11 h). Nitrocellulose 
membrane (0.2 um pore size, BioTraceTM) was laid 
onto the top of each plate for approximately 1 min. The 
membranes were transferred colony-side up to LBKan 
agar plates supplemented with arabinose (0.2%) and 
IPTG (1 mM), and incubated at 37°C for 3 h. The mas- 
ter plates were incubated for a further 3 - 5 h at 37°C, 
until the colonies reached a diameter of 1-2 mm. The 
membranes were lifted from the agar plates and placed 
on chloroform-saturated filter paper, colony-side down, 
for 1 min, after which the chloroform was allowed to 
evaporate completely. The membranes were then gently 
agitated in blocking solution (PBS plus 0.5% Tween 20 
and 5% skimmed milk powder) for 1 h at ambient tem- 
perature, washed in PBST (PBS plus 0.5% Tween 20; 4 
x 10 min) with gentle agitation and probed with the pri- 
mary antibodies (depleted antisera) in 10 ml PBST 
(1:1,250) for 16 h at 4°C with gentle agitation. The 
membranes were then washed four times in PBST and 
agitated for 1.5 h in secondary antibody solution (HRP 
conjugated to goat anti-rabbit IgG [Sigma]) (1:30,000). 
The membranes were washed four times in PBST, 
rinsed twice in PBS and washed for 10 min in PBS, 
under gentle agitation. Enhanced chemiluminescent 
(ECL) reagent was used to develop the membranes and 
the chemiluminescence was visualised by exposure of 
Roche Lumi-Film Chemiluminescent Detection Film to 
the membranes. Putative positive clones were identified 
on the master plates and each one was transferred to 
fresh LBKan agar. 



PCR verification of insert 

For verification of the presence of cloned DNA, putative 
positive colonies were used as the template source for a 
colony PCR and the T7 promoter and T7 terminator 
primers (Novagen, Notts, U.K.). Thermal cycling condi- 
tions using Taq polymerase comprised an initial dena- 
turation of 5 min at 94°C, 35 cycles of denaturation at 
94°C for 30 s, annealing at 55°C for 30 s, extension at 
72°C for 30 s kb" 1 product, followed by a final extension 
at 72°C for 7 min. 

Secondary screening 

Putative positive colonies were cultured overnight in 
BHI Kan (1 ml), at 37°C, without shaking. The cells 
were harvested by centrifugation at 16,000 g. The super- 
natant was decanted and the cells resuspended in 20 1 
BHI Kan. Each suspension was spotted in triplicate (1 
ul) onto duplicate nitrocellulose membranes and placed 
on a BHI Kan agar plate. The plates and membranes 
were incubated for 3 h at 37°C, the membranes removed 
and one of the duplicate membranes overlaid onto a LB 
Kan agar plate supplemented with 0.2% arabinose and 1 
mM IPTG while the other membrane was placed onto a 
LB Kan agar plate. These were incubated for 3 h at 37° 
C. The membranes were removed from the plates, and 
placed on chloroform- saturated filter paper for 1 min. 
Once dry, 1 ul of the lysogen-specific antiserum was 
spotted onto the bottom of the membrane, as a positive 
control. Antibody reactivity was determined as described 
above for primary screening. 

DNA sequencing 

Plasmid DNA was sequenced by GATC Biotech (Kon- 
stanz, Germany), using the T7 promoter and terminator 
primers. Sequences were translated using ExPASY's 
Translate tool http://www.expasy.ch/tools/#proteome. 
The sequences were aligned to the annotated <I>24 B gen- 
ome [GenBank:HM_208303] and CDS in-frame with the 
expression vector were documented. 

qPCR 

Induction of MC1061(<t>24 B ) cultures was performed as 
described above. A 1 ml sample was taken before addi- 
tion of norfloxacin to the cultures, and further 1 ml ali- 
quots removed at 10-15 min intervals throughout the 60 
min recovery time. RNA was immediately extracted and 
DNAse treated with TURBO™ DNAse (Ambion, TX, 
USA) according to the manufacturer's instructions. 
Absence of DNA was verified by qPCR. Each RNA sam- 
ple (300 ng) was reverse transcribed using random hex- 
amer oligonucleotides (Bioline, London, UK). Specific 
primers were designed to amplify an approximately 100 
bp region of each gene in the study (Additional file 1: 
Table SI). qPCR was performed using a StepOnePlus™ 
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Real-Time PCR System (Applied Biosystems); each reac- 
tion consisted of 1 ul of cDNA, 1 x SensiMixP/i/s SYBR 
(Quantace, London, U.K.), 200 nM of specific primers in 
a 25 ul reaction. The amplification cycling conditions 
were: initial denaturation at 95°C for 10 min; 39 cycles 
of denaturation at 95°C for 10 s; annealing at 60°C for 
30 s; extension at 72°C for 5 s. A melting curve analysis 
was performed for each amplification reaction, with a 
temperature gradient of 0.1°C from 55°C to 95°C. No- 
template controls and a calibration curve, consisting of 
6 dilutions of the PCR amplicon of each gene cloned 
into PCR-Blunt vector (Invitrogen, Paisley, U.K.) line- 
arised with Nco I (NEB, Herts, U.K.), were included in 
every experiment (Additional file 1: Table SI). Statistical 
analysis was performed using a one-way ANOVA com- 
paring gene copy numbers at different time points in 
each experiment to test the hypothesis that there is no 
variation in gene copy number during the recovery per- 
iod. A post-hoc Dunnett's test was employed, using the 
sample corresponding to the lysogen culture (-60) as the 
reference group, to assess whether or not time points 
differed from the reference. P < 0.05 values were consid- 
ered to be statistically significant. 

Protein extraction for 2D-PAGE 

Cultures of MC1061 and MC1061(O24B) were incubated 
for 6 h at 37°C. Cells were harvested and pellets washed in 
1 ml of wash solution (10 mM Tris-HCl, pH 8.0; 1.5 mM 
KH 2 P0 4 ; 68 mM NaCl; 9 mM NaH2P04). Cells were 
resuspended in 1 ml of resuspension buffer (10 mM Tris- 
HCl, pH 8.0; 1.5 mM MgC12; 10 mM KC1; 0.5 mM DTT; 
0.1% SDS; 20 ul of protease inhibitor [Roche Complete- 
Mini EDTA Free protease inhibitor cocktail tablets]) and 
each sample was sonicated for 5 x 10 s. DNase was added 
(5 ug ml" 1 ) and samples were incubated for 1 h at 37°C. 
Samples were centrifuged for 1 h at 12,000 g, the superna- 
tant recovered and protein concentration determined 
using the Bradford Assay. Aliquots (110 ug protein) of the 
sample were taken and precipitated in 10% TCA in acet- 
one containing 20 mM DTT for 45 minutes at -20°C. Pel- 
lets were washed twice in ether. 

2D-PAGE 

Isoelectric focussing was carried out on 18 cm IPG 
strips (pH 4-7,3-5.6 and 5.3-6.5;GE Healthcare), at 3,500 
V for 7 h. Proteins were separated in the second dimen- 
sion on 1.5 mm 4% stacking/15% resolving SDS-PAGE 
gels, for 6.5 h at 20 W per gel (up to maximum of 180 
W). Proteins were silver stained [50]. 

In-gel digestion of protein samples 

This was carried out according to the protocol described 
by Courchesne & Patterson [51] with the following 
modifications: protein spots were excised from the gel 



and destained with 50 ul of destaining solution (30 mM 
potassium ferricyanide, 100 mM sodium thiosulphate) 
until the silver stain disappeared; protein digestion pro- 
ceeded in 25 mM ammonium carbonate/trypsin (5 ng 
ul" 1 ) at 37°C for 16 h. 

Matrix-assisted laser desorption/ionisation-time-of-flight 
(MALDI-TOF) mass spectrometry 

Trypsin-digested protein samples were added to an 
alpha-cyano 4-hydroxycinnamic acid matrix (LaserBio- 
Labs, France) at a concentration of 10 mg ml" in 50% 
ethanol: 50% acetonitrile: 0.1% TFA. Samples were ana- 
lysed by MALDI-TOF on an ABI Voyager DE Pro 
(MALDI-TOF). The mass spectra generated were pro- 
cessed using Data Explorer to clean the spectra and iso- 
late monoisotopic peaks (all Applied Biosystems). The 
Mascot Peptide Mass Fingerprint Database was used to 
search for homologues. 

Additional material 



Additional file 1: Table SI PCR amplification primers used in this study. 
A compilation of all of the amplification primers used in this study along 
with amplification efficiency information. 

Additional file 2: Table S2. Significance of Dunnett's test results for 
gene expression data in Figure 3: Results of the Dunnett's test to 
determine significance of gene expression profile differences before and 
after prophage induction. 
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